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Abstract 

We use the quasistatic approach to analyze the criterion of ferromagnetism for 
two-dimentional (2D) systems with the Fermi level near Van Hove singularities (VHS) of 
the electronic spectrum. It is shown that the spectrum of spin excitations (paramagnons) is 
positively defined when the interaction I between electrons and paramagnons, which corre- 
sponds to the Hubbard on-site repulsion U, is sufficiently large. The critical interactions J c 
and U c remain finite at Van Hove filling and exceed considerably their values obtained from 
the Stoner criterion due to incommensurate spin fluctuations, which are important near the 
ferromagnetic quantum phase transition. Combining the quasistatic approximation and the 
equation of motion method for the Green function we obtain the results for the electronic 
self-energy to first order in the inverse number of spin components. 
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1 Introduction 



Two-dimensional and layered systems with strong electronic correlations have been attracting 
attention of researchers for more than two decades. In such systems, in particular, the 
Cu0 2 planes of the La^zSr^CuCU and YBa 2 Cu307_,5 compounds the low dimensionality 
favors appearence of the high-T c singlet superconductivity [lj. On the other hand, a number 
of low- dimensional triplet type superconducting compounds (Sr 2 Ru0 4 , UGe 2 ) have been 
discovered recently. For example, the compound UGe2 is ferromagnetic at low temperatures 
and pressure. It is expected that the ferromagnetic fluctuations may play a crucial role in the 
properties of the paramagnetic compound Sr 2 Ru04 [2]. This is confirmed by the properties 
of the La-doped compound La^Sr^zRuO^ which is near the ferromagnetic transition at 
x = 0.27 [3J, and also by the properties of the isoelectronic compound Ca 2 Ru0 4 [I], which 
becomes ferromagnetic under pressure. The above-mentioned compounds have Van Hove 
singularity (VHS) of the density of states near the Fermi level, the description of magnetic 
fluctuations in the presence of these singularities is an important problem, which may serve as 
a basis for deeper insight into magnetic and superconducting properties of these compounds. 

According to the Stoner theory, the large value of the density of states at the 
Fermi level, occuring due to VHS, leads to a possibility of ferromagnetically ordered ground 
state. The Stoner theory, however, neglects fluctuations of magnetic order parameter, which 
may significally change the necessary conditions for ferromagnetism. It is well known that 
for systems with nonsingular density of states the Stoner theory does not explain correctly 
magnetic and thermodynamic properties, e.g. it predicts larger temperatures of magnetic 
phase transition in comparison with the experimental data. Starting from the papers by 
Murata, Doniach [5] and Dzyaloshinskii [6], the paramagnon theory has been developed to 
account for the effect of spin fluctuations and describe the properties of weak- and nearly 
ferromagnetic materials. Later on this theory was systematically formulated by Moriya [TJ. 
This theory considers magnetic excitations (paramagnons) , which give the dominant con- 
tribution to thermodynamic properties of weak ferro- and antiferromagnets. Properties of 
paramagnons in the vicinity of quantum phase transitions were investigated later within the 
renormalization group method by Hertz in the framework of the 4 -model [8], the corre- 
sponding results were generalized by Millis to finite temperatures [9]. The applicability of 
the Hertz-Moriya-Millis theory in two- and three dimensions has been, however, recently 
questioned, because the lowest order corrections to the spin susceptibility were shown to be 
non-analytical with respect to momentum and therefore dramatically change the spectrum 
of magnetic excitations |10j . 

In the presence of VHS in the vicinity of the Fermi level the Hertz-Moriya-Millis 
theory itself, however, is not applicable to describe the ferromagnetic ground state. Indeed, 
due to the influence of VHS the noninteracting magnetic susceptibility Xq ma y n °t have a 
maximum at the point q = [UJ . Since the spectrum of paramagnons in the Hertz-Moriya- 
Millis theory is determined by the inverse spin susceptibility (Xg) _1 > this spectrum is not 
positively defined in the presence of VHS and the ground state with the spin density wave 
at the wavevector Q ^ becomes more preferrable in this theory. Due to the presence of 
the electron-electron interaction, however, the paramagnon subsystem is not independent, 
but it is interacting with the electronic degrres of freedom, which can change the position 
of the spin susceptibility maximum in momentum space. Since the Hertz-Moriya-Millis 
theory neglects the renormalization of the momentum dependence of the susceptibility by 
the electron-paramagnon interaction, it cannot therefore describe correctly ground state 
properties in the presence of VHS. Furthermore, for the position of the Fermi level in the 
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vicinity of the VHS the paramagnon interaction energy in the Hertz-Moriya-Millis theory 
is negative, which corresponds to attraction between paramagnons with the possibility of 
bound states formation. The magnitude of this interaction may be also changed essentially 
by higher-order terms in the electron-paramagnon interaction /, which can result in the 
repulsive paramagnon interaction. 

Therefore, to investigate the possibility of the ferromagnetic ordering in the pres- 
ence of VHS it is necessary to take into account the effect of the electron-paramagnon 
interaction on the spin susceptibility more accurately than in the existing spin-fluctuation 
theories. An attempt to go beyond these theories was made by Hertz and Klenin [12], who 
considered the sum of infinite series of diagrams for irreducible susceptibility. However, the 
renormalization of the momentum dependence of the susceptibility, which becomes crucial 
in the presence of VHS, was not considered in their study. The aim of the present paper is 
to formulate a consistent spin-fluctuation theory in the presence of VHS and to investigate 
their influence on the possibility of ferromagnetism. 

To this end, we consider the spin-fermion model, which treats interacting elec- 
tronic and spin degrees of freedom. This model was introduced in early 1990's to treat 
antiferromagnetic fluctuations in high-T c superconductors [13] (see also Refs. [Ill IS]), but 
can be generalized to the case of ferromagnetic fluctuations [IB] . Since previous attempts to 
derive this model from the Hubbard model met the problem of the so-called Fritz ambiguity, 
we use alternative way to derive this model and obtain the relation between the parame- 
ters of the Hubbard and spin-fermion models. Having this relation, the phase diagram of 
these models can be investigated. It is argued that for large enough electron-paramagnon 
interaction I c (and the corresponding Hubbard on-site Coulomb repulsion U c ) the spectrum 
of paramagnons becomes positively defined, which corresponds to the possibility of ferro- 
magnetism at / > I c . The quantities I c {U c ) exceed substantially the corresponding values, 
obtained from the Stoner criterion. The latter is in an agreement with the results of the 
functional group approach [17] and the numerical investigations for the Hubbard model. The 
critical values of the interaction I C (U C ) are finite at the Van Hove filling. 

The plan of the paper is the following. In Section 2 we introduce the spin- 
fermion model and derive the "generalized" Stoner criterion for ferromagnetism. In Sections 
3.1 and 3.2 the application of this criterion for investigation of the phase diagram and 
magnetic properties of the model is presented. The additional criterion of ferromagnetism, 
requiring the positivity of the energy of the paramagnon interaction is considered in Section 
3.3 and showed to be less important than the generalized Stoner criterion. In Section 3.4 
we combine the quasistatic approach with the equation of motion method for the electronic 
Green function which allows us to investigate electronic properties of the system. The 
obtained results are discussed in Section 4. In Appendix A we discuss the derivation of the 
spin-fermion model from the Hubbard model. The diagram technique for the spin-fermion 
model is presented in Appendix B. The formulae for the second derivative of magnetic 
susceptibility with respect to the momenta are presented in Appendix C. 
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2 Magnetic Properties of Two- Dimensional Systems 



2.1 The Hubbard model and an applicability of the Stoner crite- 
rion of ferromagnetism 

To investigate the effect of VHS on magnetic properties of 2D systems we consider the 
Hubbard model 

U = /J ^c^ a c ka + U 22 n it n ii- (!) 

ktr i 

with the electronic dispersion 

£k = — 2t(cos k x + cos k y ) + 4t'(cos k x cos k y + 1) — \i. (2) 

where t and t' are the nearest- and next-nearest neighbor hopping parameters, U is the 
Hubbard on-site repulsion, n ia = c[ a c ia . 

Let us first consider the noninteracting case. The corresponding plots of the non- 
interacting density of states for different t' are shown in Fig. [TJ For t' < 0.5t the density of 
states has logarithmic singularity at the energy At — 8t' measured from the bottom of the 
band. For t' = 0.5t there are lines of VHS along the k x = and k y = directions, and 
the density of states has stronger divergence p(e) oc l/e 1 ^ 2 at the bottom of the band (the 
so-called flat band case [T7J HE]), which is analogous to the giant VHS in 3D systems [19|. 
The dynamical susceptibility of the noninteracting electronic gas with the dispersion (J2J) is 

Xq = -Jj^GkGk+q (3) 

k 

where Qk = i}v n — ^k) -1 is the bare electronic Green function, k = (k; iu n ) and q = (q; iu> n ), 
v n = (2n + 1)ttT and uo n = 2nnT are the Fermi and Bose Matsubara frequencies. In the 
presence of VHS the uniform susceptibility is logarithmically divergent near Van- Hove band 
fillings with /i = 0, 

Xo — m 7T1 7- ( 4 ) 

2nt max(|/i|, T) 

According to the Stoner criterion, this leads to a possibility of ferromagnetism for U > 
jystoncr _ 1/^-0 However, the momentum dependence X( q -o) a ^ ^ ow T <^t and \i 7^ has its 
maximum at q 7^ (see Ref. [11] and Figs. I4f5l with A = below). In the random phase 
approximation (RPA) 

*T - T% q w 

this results in an instability of ferromagnetic ground state, since xf^o) is divergent at q 7^ 
even for U < U^ toncT , which implies the spin-density wave instability in this approach. This 
shortcoming holds also in the Hertz- Moriya-Millis theory, not considering the renormalization 
of momentum dependence of the electron-electron interaction and implies the impossibility 
of ferromagnetism near van Hove fillings in existing spin-fluctuation approaches. 

2.2 The spin-fermion model 

For the purpose of investigation of magnetic properties beyond RPA-type approaches we 
consider the spin-fermion model 



ZAV,V + ] = f D[S ]C ,c + ]exp(-(3S), 
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Figure 1: Electronic noninteracting density of states p(e) at a) t' — 0.45£, b) t' = 0.50i 



5 = Yl C t(- iz/ n + £ k)c fc(7 + S S + S <? S -9 - Yl (Vtck* + 4aVka) , (6) 

ka q ka 

where the Grassman fields Ck a , c~l a correspond to electronic degrees of freedom, the field 
S q describes the collective spin degrees of freedom, corresponding to paramagnons, Xg is 
the "bare" susceptibility of the paramagnon subsystem, s,j = J2kaa' ^ta^^' c k+q^' is the 
spin operator of itinerant electrons, a are the Pauli matrices, I is the interaction of the 
paramagnon and electron subsystems, i]ka,Vka are the fermionic sources fields, j3 = 1/T is 
the inverse temperature, and N is the number of sites. 

Although the model was previously derived from the Hubbard model via the 
Hubbard-Stratonovich transformation, this way of derivation meets the problem of the so- 
called Fritz ambiguity [3U1 ED E21 E21 EH ES] and the resulting action is not uniquily defined. 
In this paper we use a different way to establish connection between the spin-fermion and 
Hubbard model, and show that the former describes magnetic degrees of freedom of the 
latter, see Appendix A. The equivalence with the Hubbard model requires 



U, Xq 1 = (xtr 1 + U 2 H, 



q- 



(7) 



where xf is the exact magnetic susceptibility of the Hubbard model, which is calculated 



self-consistently in the framework of spin-fermion model according to (x[ 



H~ l -U,H q 



is the irreducible electronic polarization operator. In the most part of the paper we consider, 
however, the spin-fermion model with arbitrary Xq and /. 

To investigate the magnetic properties of the model it is convinient to rewrite 
the action in terms of the bosonic fields only [12]. Integrating out electronic degrees of 
freedom and expanding the resulting functional in the series in / we obtain 



D[S) exp 



kk', 



:\. j 

Sint = T } -Tr 



n=2 



N 



Qk^k-l 



ka,k'a' 



where 



^ka,k'a' [S] — {{iiO n — e-k)5kk>5aa> ~ {I / ' N)B aa iSk-k'} ka k'a' 



(9) 
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is the Green function of electrons, propagating in the presence of an external field S, S- mt 
corresponds to paramagnon interaction, the exponent n corresponds the power of matrix in 
the square brackets taken with respect to indices ka, k'a'. Coefficients of the expansion of 
the interaction iS int in powers of the spin fields S 3 k _ k , (j = x, y, z) determine the vertices of 
paramagnon interaction 



TF 



i=l 



a 



30 



(10) 



(V g corresponds to all possible permutations of momenta gj, Tr^ is the trace with respect 
to spin variables). Taking into account the fourth order vertex only leads to an effective 4 
model [H [9] with the paramagnon dispersion (x^ PA ) _1 an d the interaction r (gi, q2, <?3, #4)- 
Below we consider, however, the effect of infinite sequence of vertices ( TTUj) on the magnetic 
properties. 

Diagram technique for the model ([S]) is considered in Appendix B and contains 
the following elements: 



• the propagators of the longitudinal and transverse paramagnons 7Z q ; 

• the electronic Green functions Qk, connected in loops (1TU1) ; 

• the electron-paramagnon interaction vertices la. 

These rules are used below to calculate physically observable quantities. 



2.3 Spin susceptibility 

To investigate magnetic properties of the model §6§ we consider magnetic susceptibility, 
defined as the causal double-time Green function 



1 



47V 

and the paramagnon propagator 



4' = 7X7 / (T r <(r) S=q (0)) (11) 







K* = Jff d ^ nT <TrSj(r)Sl q (0)>. (12) 

Here s^(r) and S^(t) are the operators in the Heisenberg representation, corresponding to 
the and fields, T T denotes the chronological product with respect to the imaginary 
time. To obtain the result for the susceptibility xf, one can differentiate Eq. (jSJ) with respect 
to rj. Classifying corresponding contributions to the one particle reducible and irreducible, 
we obtain (see Appendix A) 

H, 



Xt^^j^ + Xt (13) 

where xf is the non-paramagnon contribution to the electronic susceptibility (xf = UH^ 
for the Hubbard model). The result (|T3|) reduces for H q = x q to the RPA result for the spin 
susceptibility. 

The general expression for the paramagnon propagator is 

K- 1 = Xq 1 ~ S ? , (14) 
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Figure 2: (a) Series of diagrams for the polarization operator of electronic subsystem; the dashed circle 
represents the sum of all one-loop diagrams, as in Fig. 2b. (b) The example of the one-loop diagram, double wavy 
lines correspond to the longitudinal or the transverse propagator, denoted by dashed and wavy lines, respectively 

where S g is the paramagnon self-energy. It may be established diagrammatically that this 
self-energy can be expressed through the polarization operator H q according to E q = I 2 H q 
(this relation is in fact a consequence of the form of the interaction term of the model (jBD). 

The polarization operator H q can be represented as a set of diagrams, which con- 
tain loops of electronic Green functions (the one-loop vertex functions, defined in Appendix 
B), connected by two or more paramagnon lines (see Fig. 2). In the present paper we neglect 
contributions to H q , which contain more than one electronic loop (second and next terms 
in Fig. [2^), since we expect that their contribution to the magnetic susceptibility is small. 
The contribution of the diagrams with only one electronic loop (Fig. 2b) corresponds the 
one-loop 2-point vertex 

H q = (njS]> =: -^Lrf_ loop (g, -q), (15) 
njS] = ~ J2 Tl ° [G fcl , fc2 [S]a z G fc2+9ifcl+9 [SK] . 

fci/C2 

According to the Mermin- Wagner theorem [21], the long range order in 2D sys- 
tems is possible at T = only. The susceptibility (fT51) in the limit T —>■ corresponds to 
the susceptibilities of the ordered state, averaged over directions: 

hm X e q l = l(x z q z + xt) C T=o- (16) 

The necessary condition of the ferromagnetic ground state is 

I 2 X H q=0 = 1 at T -> (17) 

(generalized Stoner criterion), and the positivity of the spectrum of excitations of the static 
paramagnons u; q = 7^(q. ) at T ^ 0, which is fulfilled if the product X(q;o)-#(q,o) is maximal 
at q = 0. This criterion is violated in RPA, where u; q is not positively defined in the T — > 
limit in the lowest (second) order of perturbation theory in I (see Section 2.1). Moreover, 
the corresponding energy of the paramagnon interaction To(0, 0,0,0) = — 12J 4 p"(/i) is also 
negative near VHS. To investigate the possibility of ferromagnetism we consider below the 
results for the polarization operator H q beyond RPA. 

2.4 The static and the quantum contributions. The static approx- 
imation 

In the assumption of ferromagnetically ordered ground state considered expressions can be 
simplified in the T — > limit analogously to the earlier investigated case of antiferromagnetic 
order [H [151 US] . 
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Figure 3: The diagrams for the polarization operator H q in the lowerest order in /, the notations are the same 
as in Fig. 2 

Let us consider for example the lowerest (second) order result of perturbation 
theory in I for the polarization operator susceptibility H q . There are two topologically 
different diagrams in this order, see Fig. [3J The corresponding analytical expression reads 
(cf. Appendix B) 

6T 2 I 2 T 2 I 2 

^<? 2) = J^T~ E ^kSk+qSk+g'T^q' + ~^T2~ E 'SkQk+qQk+q'Qk+q+q'Tlq' , (18) 

kq' kq' 

We assume that the paramagnon propagator lZ q at small momenta and frequencies for the 
Fermi level lying near VHS has a form (see, e.g., Ref. [22] ) 

^ = q 2 + ^ 2 + B\u n \/ max(tg + g_, \u n \) + (19) 

where A, B > are some constants, £ is the correlation length of spin fluctuations, 7^q i(Jn is 
the contribution, which is regular at q — > 0, £ — > oo and neglected below. The result (fl9l) 
can be obtained from the general formula (THI) under the assumption that the expansion 
of the polarization operator H q in the momentum and frequency has the same form as one 
for the noninteracting susceptibility Xqi the effect of the interaction is in the change of the 
parameters A and B only. While the condition B > being a consequence of the analytical 
properties of the susceptibility as a function of frequency, the positivity of the constant A is 
analyzed below in the section 3.1. 

It follows from the Eq. ( fT9l) that the dominant contribution to the momenta- and 
frequency sum in Eq. ( |T8|) comes from the region with |q| ~ and the zeroth Matsubara 
frequency. If the condition 

{t/Tf' 2 <^i (20) 

is fullfilled the contribution of the terms with nonzero Matsubara frequency in the sum 
over uj' in the Eq. (fl9|) can be neglected, together with the q-dependence of the electronic 
Green function. The condition (T2"U|) is certainly fulfilled in 2D case at finite T above the 
ordered ground state, since of the correlation length is exponentially large at small T (cf., 
for example, Refs. [231 121] )• This differs present theory from the earlier considered 3D 
case [12J, where abovediscussed approximations give only qualitatively, but not quantitavely 
correct description of magnetic properties. 

Within these approximations the sum (|18|) takes a form 

H ? ] = ^? E (~ 6 ^ + g&U E -f- ( 21 ) 

N k q' "V 

Analogous approximations are also applicable to higher order diagrams. The results coincide 
with those for the model with the action 



Z^oc [r],r] + } = J d 3 S exp 



3I 2 S 2 
2A 2 



P ^ Vkfon - £k - ifrS) v 



(22) 



which contains only one uniform static mode S with the propagator 



l(A//) a = lri 



(23) 



The Green functions of Bose and Fermi fields of the model (1221) can be expressed 
in closed analytical form, cf . [T5], [16] . The static approximation and the neglection of q- 
dependence of electronic Green function leads to zero momentum- and frequency transfer 
along the paramagnon lines in all diagrams; hence, Gfcfc/[S] becomes diagonal with respect 
to momenta and frequency: 



<GW[!3] -> G k (S)5 kk >, G k (S) = (iu n -s k - IaS) 



-i 



Therefore we obtain 
where 



n,(S) 



H 9 = (iyS))^ 
T 



2N 



^Tr ff [G fc (S)<7*G fc+ ,(S)<7*] , 



(24) 
(25) 
(26) 



the subscript £ — > oo corresponds to averaging with the functional (J22J). Calculation of the 
average in the Eq. fl25|) leads to 



d 3 S 



~IL\\(q\S signS*) + ~n ± (q|Ssign^ 



exp ( 


3I 2 S 2 \ 




2A 2 J 



(27) 



where 



n||,±(q|5) 



N 



f(e k -IS)-f(e k+cl TlS) 



Similarly we find the electronic Green function: 

5 2 Z si [r),r) + 



IS±IS 



Gi 



(G*(S)> 



f28l 



The explicit expression for Gk was obtained earlier in the paper [16]. The corresponding 
spectral function 



A(k,u) = —ImG k (u) = JL 



v cj-e k ) exp 



3(u; - e k ) 
2A 2 



(29) 



has two-peak (non-quasiparticle) structure at the Fermi surface, which destroys the quasi- 
particle picture due to strong ferromagnetic fluctuations (see Fig. [10] below). As it was 
discussed in the paper [23], the corresponding violation of Fermi-liquid behavior corresponds 
to a quasi-splitting of the Fermi surface at low temperatures, which is related to the change 
of the the electronic spectrum in the vicinity of the magnetically ordered ground state. 



3 The Results for the Phase Diagram, Paramagnon 
Vertices and Spectral Functions 

3.1 The irreducible susceptibility 

To investigate the possibility of the ferromagnetic order in the ground state we consider the 
momentum dependence of the irreducible static susceptibility (polarization operator) i?( q ,o) 
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Figure 4: Momentum dependence of the static irreducible electronic susceptibility iJ q in the first quarter of 
the Brilluin zone for t' Jt = 0.45, n = 0.583 > «vh = 0.466 and different values of A 



A/t=0.0 A/t=0.1 A/t=0.2 




Figure 5: The same as in Fig. |U for n = 0.338 < n V H 

at T — > 0. The plots of -ff( q ,o), calculated according to the Eq. ([25]) in a quarter of the Brilluin 
zone for t' = 0.45t and different values of n, A are presented in Figs. HI [5] The chemical 
potential /i is adjusted to keep the number of electrons n = (2/N) f_ de A(k, e) equal 
to the noninteracting value with increasing A, the spectral function A(k, e) is determined 
by the Eq. (j25). 

For A = (which corresponds to / = 0) the global maximum of the static 
susceptibility is located at the point q ^ 0, so that for both the positions of the Fermi level 
above and below the VHS the condition of the positivity of paramagnon excitation spectrum 
is violated. With increasing A the static irreducible susceptibility has a maximum at q = 0: 
for n < nvH the local maximum at q = becomes the global one and for n > nyH local 
maximum shifts to the point q = 0. Therefore, in both cases the interaction of electrons 
with paramagnons leads to the global maximum of polarization operator of the electronic 
subsystem at q = 0, which corresponds to a possibility of the ferromagnetically ordered 
ground state. 

To determine the critical value of A (the minimal value of A, at which the fer- 
romagnetic ground state is possible) we consider the second derivative of the irreducible 
susceptibility with respect to q x (or q y ) d^H^ = d 2 H cl /dql at the point q = (cf. Ap- 
pendix C). The results of calculation of this are presented in Fig. [6j Changing of sign of 
the second derivative determines the critical values A c ~ t, which depend on the electronic 
concentration. 
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Figure 6: The derivative d% if q= o vs A for: n — 0.533 — solid line, n — 0.583 — dashed line, n — 0.626 — 
dash-dotted line 





Figure 7: Phase diagram in n — U plane: (a) t'/t = 0.45; (b) t'/t = 0.50. Solid line — the result of the 
present paper; dashed line — the results of functional renormalization approach methods [17j; dash-dotted line 
— the result of the mean field theory 
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3.2 The phase diagram of the Hubbard model 

To obtain the phase diagram it is necessary to establish the relation between the critical 
value of parameter A c and the corresponding spin-fermion interaction J c or the Hubbard 
interaction U c . These relations can be found using the generalized Stoner criterion (I17p . 
Below we consider the Hubbard model, for which we find U c = 1/H (A C ), where H is the 
uniform static polarizaton operator, 

4+8t' 

Ho = ^r [ P(e)v( £ -^)de, (30) 



2A J rv ^ \ A 

-4+8t' 

(p(x) = (3x 2 + 2) exp(-3x 2 /2)/ v / 6^. 

The resulting phase diagram in the n — U plane for t'/t = 0.45 and n > wvh 
is presented in Fig. [7^ (the case n < n V H is not considered in the following, because 
the incommensurate spin fluctuations with the wavevector Q far from the point q = 
are not taken into account in the present approach). The critical values U c for n > tt-vh 
are larger than the corresponding mean-field values . The results for n far from wvh 
are in qualitative agreement with the functional renormalization group (fRG) results for 
the Hubbard model [17]. However, in contrast to the fRG results the critical value U c is 
nonzero for n = nyH, which is related to the breakdown of the quasiparticle picture of the 
electronic spectrum owing to the ferromagnetic fluctuations. Thus, the non-Fermi-liquid 
properties of the electronic spectrum in 2D systems become important for the criterion of 
the ferromagnetism in the vicinity of the VHS. 

The phase diagram for the flat band case t'/t = 0.50 is presented in Fig. [7b. In 
this case the critical values U c are also larger than the mean field values , and there is 
qualitative agreement between the fRG and obtained here results for all values of electronic 
concentration n 1171. 



3.3 The paramagnon interaction vertex 

Let us consider the vertex T 333 ' 3 1 (qx, q 2 , <?3, which determines the energy of the param- 
agnon interaction and defined by the Eq. (|47j) of Appendix B. Below we are interested 
in j ^ j'-the component of this vertex. This component can be expressed through the 
corresponding 2-particle irreducible vertex 3 (qx, q2, 93, q<i), which is defined such that it 
cannot be split into two connected parts by removing two paramagnons lines. In general the 
relation between these two vertices is determined by the parquet set of diagrams. 

To simplify the relation between the vertices we generalize the model ([6]) intro- 
ducing the M-component spin field S (cf. Refs. [25l 126]). j = 1, . . . , M. Due to internal 
sums over spin indices, the dominant contribution into the full vertex in the limit M — > oo 
comes from the ladder set of diagrams with equal spin indices of internal paramagnon Green 
functions (see Fig. [8]). 

This allows to obtain a simple relation between the full and irreducible vertex, 
V 333 3 and 3 to first order in 1/M 

T 333 ' 3 '(qx, q 2 , g 3 , Qa) = T 3 . 33 ' 3 ' (q x , q 2 , q 3 , qi) 

~ giv rJjy ' y '(^'?2,p,p0^p^p' r C y ' yy (-p ) -p' ) 93,g4)5 (?1 +< ?2 ,-p-p'- (31) 
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Figure 8: The ladder diagrams for the paramagnon vertex 
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Figure 9: The irreducible vertex (t 3 /I )Ti IT vs A at different electronic concentration n and t'\ (a) t' — 0.45i: 
n = 0.533 — solid line, n = 0.583 — dashed line, n = 0.626 — dott-dashed line; (b) t' = 0.5t: n = 0.316 — 
solid line, n = 0.415 — dashed line, n = 0.486 — dott-dashed line 



At large values of the correlation length £ the dominant contribution to the sum (13"T1) comes 
from the small momenta p,p'. Supposing that the external momenta (?2, ?3, ?4 are also 
small and Y\" 1 (qi, q2,Q3, 94) is constant in the vicinity of qi = q2 = q3 = q4 = (below we 
denote this constant as r; rr ) we obtain in the framework of the quasistatic approach 



where 



Y n3 ' 3 '{q 1 ,q 2 , q^qi) 



p 



;^i,0"w2,0"^3,0"a;4,0) 



l + r irr MC qi+q2 /2 

TA 2 . |q| + v/q 2 + 4£" 2 



4vr|q|Vq 2 T4e z 2 



In 



2£ 



-1 



(32) 



(33) 



Note, that in the limit £ — > 00 the quantity £ q = TA 2 /(4:ir\q\ 2 ) ln(|q|£) diverges at |q| — > 0, 
which results in vanishing of the full vertex of paramagnon interaction (1321) at q x = q 2 = 
q 3 = q 4 = in accordance with the assumption of existing of long-range order in the ground 
state q 

One can see from the Eq. ( 1321) . that the sign of Fi rr at small external momenta 
is crucial for the possibility of the FM gflound state. If r irr < 0, the expression (152]) shows 
that paramagnons can form bound states, and hence the ferromagnetism is not possible. 



1 Actually this vanishing of the vertex generalizes the Adler principle [57] for the electronic system in- 
teracting with paramagnons; this is similar to vanishing of the magnon- interaction vertex at qi = in 
Heisenberg magnets which was discussed earlier (28] . 
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Let us consider the sign of r irr within the quasistatic approach depending on the magnitude 
of the electron-paramagnon interaction. In general, the vertex r irr can be represented dia- 
grammatically as a sum of the contributions with the one-loop vertices r 1 _i oop , connected 
by three or more paramagnon lines. Using the same approximation as for the one-particle 
polarization operator r irr = ri_i oop , i. e. neglecting contributions to the irreducible vertex 
with the number of loops n > 1 we obtain in the M = 3 case 

T irr = ^Tr CT (^G k (S)a z G k (S)a*G k (S)a z G k (S)a x \ (34) 

= h III, p{£)9 i 1 ^) d£ ' g{x) = i^ (2 + 3x2 ~ 9x4) exp ( _3x2 / 2 ) • 

The results of the numerical calculation of the irreducible vertex in the T — > limit using the 
Eq. (1341) are presented in Fig. 9. One can see that the irreducible vertex r irr changes its sign 
at the values of A c , which are much smaller than the values of A c , obtained in Section 3.1, 
and hence they do not result in additional limitations on the possibility of the ferromagnetic 
ground state. 

For Tin- > and small momenta (Zi , • • • , ?4 one can neglect the unity in the 
denominator of (1321) in the limit M — ► oo to obtain 

T 3 ' 33 ' 3 (q u q 2 , g 3 , = -tjj. ^1,0^2,0^3,0^4,0- (35) 

iW Sqi+q 2 

Note, that for T- 1TT <C 1 (i.e. in regime of weak coupling and/or smooth enough density of 
states) the formula (135j) is correct in the considered limit of low temperature (£ — ► oo). For 
higher temperature the unity in the denominator (l32l) has to be retained also in the limit of 
large M, since T iTT ( «C 1. 

3.4 The electronic self-energy and spectral functions 

The first order 1/M result (|35|) for the paramagnon vertex can be used to investigate the 
influence of the paramagnon interaction on the electronic self-energy S and other properties 
of the electronic system. In the framework of the 1/M-expansion it was shown [26] that the 
electronic properties of the model ([6]) to first order in 1/M are determined by the self-energy 
Sfc and one-particle irreducible vertices of electron-paramagnon interaction 7^ 7^, 7^ . In 
the limit £ — >• 00 these quantities can be obtained from the system of the algebraic equations 

Sfc = A 2 7 fe G fc , 

lk = 1 + ( T f G fc - riGl) A 2 /3, 

7 f = A 2 (2 7fc 3 ^. + lkl z k z G 2 k + 7f ± G fe ) + a lk G k) (36) 

7f ± = _2A 2 ( 7fc 4 G 4 fc + 7fc 2 7f G 3 fc + lkl z k z± Gl) - ai lGl 

where a = MV 333 3 '(0, 0, 0, 0)(o- The results of the solution of these equations with account 
of paramagnon interaction (a = 2 from (1331) ) and without it (a = 0) are presented in Figs. 
6-8. 

As well as in the absence of the paramagnon interaction [16] . the real part of 
the self-energy has an infinite slope at the Fermi level, and imaginary part has a 5-function 
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Figure 10: The spectral density A = — -ImG vs uj = v — e^ + ii, calculated from (|36p without (a = 0, dashed 
line) and with account of (a = 2, solid line) the paramagnon interaction ; dash-dotted line — the same in the 
quasistatic approach (the formula ([29]) ), A = O.lt 
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peak. It was pointed out above that these anomalies originate from the electron-paramagnon 
interaction and violate the quasiparticle picture in the vicinity of the Fermi level. Physically 
this corresponds to appearance of the new quasi-split-Fermi-surfaces shifted from the Fermi 
level by ±A in the presence of strong ferromagnetic fluctuations [21]. Note, that existence of 
the window where the spectral weight is zero is an artefact of 1/M-expansion: the comparison 
with the results of the quasistatic approach shows (Fig. [TO]) that the spectral weight in this 
window is reduced, but is not zero. In the presence of the paramagnon interaction the 
spectrum is less coherent than without it (ImS is enhanced, A(u) is smeared), the distance 
between the peaks of the spectral density A is being reduced. Thus, account of the interaction 
is important for the description of the electronic spectrum. 

4 Conclusions 

We have considered the necessary conditions of the existence of ferromagnetism in 2D systems 
in the presence of VHS: the maximum of the polarization operator at q = 0, which ensures 
positive defmiteness of the magnetic excitation spectrum and the repulsion between the 
paramagnons (the latter guarantees the impossibility of the bound state formation). 

For the Fermi level position in the vicinity of VHS and in the absence of the 
electron-paramagnon interaction / these conditions are violated: the momentum dependence 
of spin susceptibility has its maximum at the wave vector q ^ 0, which corresponds to 
possibility of the spin density wave in the ground state and the paramagnon interaction is 
attractive. However, with increasing / the static spin susceptibility acquires a maximum at 
q = and the paramagnon interaction changes its sign. The critical value of the electron- 
paramagnon interaction I c (or the corresponding Hubbard interaction U c ) for which the 
ferromagnetic ground state is possible exceeds substantially the corresponding value 7^ IF 
(t/ c MF ), determined from the Stoner theory, and agrees with the functional renormalization 
group investigations of Hubbard model for the electronic concentration being in the vicinity 
of the Van Hove filling. The critical values I c and U c are however finite at the Van Hove 
electronic filling for t' < 0.5t, which is related to the non-quasiparticle picture of the electronic 
spectrum (quasi-splitting of the Fermi surface). 

In contrast to the Moriya theory, where the difference of the conditions of fer- 
romagnetism from the Stoner criterion originates from the quantum contributions to the 
uniform irreducible susceptibility, in the present approach it results from the change of the 
momentum dependence of the susceptibility by classical spin fluctuations. The change of the 
uniform irreducible susceptibility by the interaction is also taken into account in the present 
approach, but it is less important. The kind of the transition from the ferromagnetic to the 
paramagnetic phase is an open question. One of the possibilities is the existence of interme- 
diate state with the strong short-range order, characterized by the wave vector Q ^ 0. To 
describe magnetic properties of this phase, it is necessary to consider the excitations with 
incommensurate wave vector Q and account for quantum fluctuations, which is a problem 
for further investigations. The problem of ferromagnetism formation at n < nyH requires 
also further consideration. 

Besides considered reasons of deviations from the Stoner criterion, the effect of 
screening due to particle-particle scattering [2S] may be important in the Hubbard model. 
The agreement of the obtained results with the renormalization group results for this model 
at fillings away from the the Van Hove filling shows, however, that the change of momen- 
tum dependence of susceptibility seemingly plays a dominant role. On the other hand, the 
importance of the non-quasiparticle spectrum of excitations in the vicinity of ferromagnetic 
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instability, which has not been taken into account by previous fRG considerations, is demon- 
strated in the framework of our approach. 

The obtained results may serve as a basis for further investigations of itinerant 
systems in the vicinity of the ferromagnetic instability in the presence of VHS in the electronic 
spectrum. More detailed investigations of the vicinity of the quantum phase transition and 
the role of nonanalytic corrections in the system with presence of VHS are needed. It is also 
of interest to investigate a possibility of the triplet pairing in systems with of VHS in the 
electronic spectrum. 
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Appendix A. The Derivation of the Spin-Fermion Model 



In this Appendix we consider the derivation of the spin-fermion model starting from the 
Hubbard Hamiltonian ([T]). 

The originally proposed way of derivation of this model [HI [32] was to use the 
decomposition of the interaction term into charge and spin degrees of freedom, e.g. 



Un^n, 



U 



(37) 



(the summation over repeated spin indices is assumed) with the subsequent decoupling of 
quartic terms of fermionic operators using auxiliary scalar or vector fields. The representation 
( 1371) is, however, not unique, in particular another, SU(2) symmetric representation, 



Un^n il 



U 
1 



(cf a c ia ) 2 5 aa , - ^(c+a ua ,c i(T ,) 2 



(38) 



may be used (see, e.g. Ref [7]). 

It was previously discussed (see, e.g. Refs. [30l EU E2J, [33l [351 ES]) , that although 
the representations (1371) and (!38|) are equivalent, they lead to different effective actions (the 
so-called Fritz ambiguity). In particular, the representation (137)) leads to an action with the 
scalar auxiliary fields and therefore is not transparently SU(2) invariant. The representation 
0381) is free from this problem, but does not reproduce correctly the mean-field results for 
the Hubbard model due to an extra factor of 1/3. The reason for this difficulty is that 
the factor 1/3 in Eq. fl38l) takes into account that not only a lowest order component of 
the longitudinal interaction (c^a^,^/) 2 , but also the higher order transverse fluctuations 
which arise from the (cf a a^,c ia ') 2 part of the interaction contribute to the self-energy of the 
S z field and similar to the other quantities (see, e.g., Ref. [53]). 

Below we construct an effective action which is explicitly SU(2) invariant and 
at the same time reproduces correctly the mean-field results. Although some attempts of 
derivation of such an action were performed earlier, they were either restricted to systems 
with additional orbital degeneracy [32] or used additional auxiliary vector field to perform 
averaging over directions in Eq. (I3T1) [30, 34J. 

We start from the Hubbard model action 
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= ^(-i^n + ^)4a°ka + <Sint ~ ^ ^ka C ka + 4aVkc) 



ha 



S\vk — Jjy '^2(n q 7l- q - S Z q S Z _ q ), S q — C ko~<? av' C k+q<r'i 

q kaa' 

where we have introduced fermionic source fields i],7] + . To derive an effective action we 
introduce auxiliary field S q such that 



(39) 
(40) 



n q n- q + U eS x„ (s g s_ 



(s z q s z _ q ) 



(41) 



where R q 1 



(x 



0\-l 



U e ff,U e s = i7/(l + U\^),X q is th e bare polarization operator (see 



Eq. Q). It can be shown that the integration over the field S q leads to the original action 
fj4*0l) . Let us now consider the sum of the ladder diagrams, generated by the interaction 
U e ffS q St q . This summation results in the renormalizaton of the z-component of the inter- 



action itself U e ss z s 2 _ 



Us z s z _ , the electron-paramagnon vertex, U e $s z Si 



the renormalization of the paramagnon propagator R q — > Xqi where 



USqSi q and 



Xn 



(xiir'-ua-ux"). 



(42) 

where the index z-irr means 



Introducing the two-particle irreducible interaction (s q s z _ q ) z -i TT , 
that the ladder diagrams generated by this interaction have to be excluded, the result (14TI) 
can be therefore rewritten in the form 



= jf E {x^s z q s% + R-'is^s^ + s»s%) + u cS ( s *si q + s y q si q ) . 

q 

+Us z q S% + j [n q n_ q - (s z q s%) z ^ r ] + U eSX ° q (s*s*_ q + S y_ q ) 



(43) 



Note, that in addition to the spin-fermion interation the action (1431) contains the Hubbard 
interaction irreducible in the s^-spin channel. Performing the same steps in x and y spin 
directions we obtain 

Sint = Yl {^^ S </ S -? + us q s -q + J [n 9 n- q - (s 9 s_ 9 ) ph _ irr ] 1 (44) 
q 

where the subscript ph-irr denotes that the ladder diagrams in longitudinal or transversal 
particle-hole channel generated by the interaction (s 9 s_ g ) are excluded. Consideration of 
higher order in U contributions to the ladder-type processes leads to the same action fj44l 
with the exact polarization operator H q instead of the bare one x q i n the Eq. fj42|) . 

The effective model flH|) contains both the spin-fermion (first two terms) and the 
irreducible Hubbard (third term) interaction. Due to its irreducibility the latter is expected 
to lead only to small renormalization of the parameters of the spin-fermion model and it is 
neglected in the main text. 

To obtain the magnetic susceptibility of the Hubbard model (|39|) we differentiate 
the corresponding generating functional (1391) . fj4*4"|) with respect to fermionic sources. In this 
way we obtain two contibutions: 

X e q =X e q ,sF + UH 2 q , (45) 
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where the first term originates from the one-particle reducible diagrams containing the spin- 
fermion interaction US q s_ q , 

= tJ^' (46) 

and the second term originates from the irreducible Hubbard interaction. The result fT46]) is 
used in main text, Eq. (fl3|) . 



Appendix B. The Diagram Technique for the Param- 
agnon Vertex in the Spin-Fermion Model 

To consider spin fluctuations in the spin-fermion model (E]) it is convenient to introduce the 
r -vertex of paramagnon interaction, 

r jl - Jr (gx, . . .,q r ) = -2T (j[x^ J d n . . . cir r e^ + --- +i --(T r [^(r 1 )...^;(r r )]) c , 

(47) 

gi + . . . + q r = 0, 

where q±, . . . , q r are the paramagnon momenta, ji, ■ ■ ■ ,j r are the spin indices. The product 
m=i Xq~i removes the external Green functions of the field S; for r > 2 we assume x qi — > lZ qi 
which corresponds to removing the exact propagators of the field S, the average with the 
index c denotes that the diagrams, for which the external spin operators (14?]) are connected 
with each other through the lines corresponding to the electronic Green functions (see be- 
low), should be considered. Whereas the 2- vertex describes self-energy corrections to the 
paramagnon propagator, the higher order vertices describe interaction of the renormalized 
(physical) paramagnons. 

The calculation of correlation functions of spin operators can be carried out using 
the expansion exp[— [3S- m t\ of the functional (jBJ) in the series in electron-paramagnon interac- 
tion. According to the Wick theorem, the average in (T4"7]) is expressed through all possible 
combinations of averages of pairs of operators 



S 3 q S 3 ql ) =T I drdr' e^ T ^ ,T (T T 



o 



Si(r)SC(r f ) 



-Nll q 5 q „ q ,5 jf , 



where lZ q is the paramagnon propagator. The corresponding contributions to the vertex can 
be represented by diagramms, every diagram for the r-vertex (T4T1) of the order n consisting 
of: 

• n internal {a ± ,a z ), r external (S ± ,S Z ) interaction vertices, and n solid lines, corre- 
sponding to the electronic Green functions. One electronic line enters and goes out of 
each internal vertex; 

• the longitudinal vertices (a z , S z ) are connected by dashed lines, the transverse vertices 
(o~ ± ,S ± ) are connected by wavy lines. The external S* + -vertex is connected with the 
internal cx + -vertex, the ^"-vertex with a", and 5* 2 -vertex with a z ; internal wavy lines 
connect a + - and <7~-vertices, and dashed lines cr z -vertices; 
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• The factor (l/2)TZ q corresponds to each internal dashed or wavy line with the 4- 
momentum q, the factor — to each solid line with the 4-momentum k; the sum 
of momenta, which enter each internal vertex, is equal to the sum of momenta go- 
ing out of this vertex, the summation over all independent momenta in the resulting 
expression being performed; 

• The number of internal transverse vertices is even for each loop and <r + - and <T _ -vertices 
alternate along the electronic Green function line. The type of all transverse vertices in 
diagram for which it is possible to determine it unambiguously should be determined; 

• The factor corresponding to a diagram is I n (T / N)^ n ~ r ^ 2+1 (— l) c +/ +1 2 n±+nx , where c is 
the full number of longitudinal vertices, placed between neighboring transverse a~- and 
<7 + -vertices. / is the number of electronic loops, Uk is the number of the independent 
internal transverse vertices, type of which is not determined unambiguously, n± is the 
number of wavy lines; 

The lowest-order diagram for the vertex of any order is the electronic loop (|10p . 
The vertex renormalization is performed in two ways: the internal paramagnon lines can 
appear in electronic loop or consist of several loops, connected by paramagnon lines. 

To classify diagrams it is convenient to define the one-loop vertex 
^l-ioapC? 1 ) •••>?»•) as a sum °f diagrams for the vertex containing only one electronic loop of 
electronic Green function (see Fig. 2b). Every diagram for the vertex consists of one or more 
one-loop vertices, connected by paramagnon lines. Comparing coefficients of expansions in 
the series in / of (T4T|) and (l4"8"j) . it is possible to show that 

r{ 1 :r^ P (9i,-..,^) = ^E T ^ E 

(GkiM [ S K lG fc2+<n,fc3+<n [SK 2 • • • G kr+qi+ ... +qr _ lM+qi+ ... +qr _ 1 [S}a Jr ) s , (48) 

index V q points that the sum is taken over all permutations numbers i of momenta and 
simultaneously spin indices jV The explicit expressions for one-loop 2- and 4- vertices are 
given by formulas (1271) and (1541) of main text. 



Appendix C. The Second Derivative of the Magnetic 
Susceptibility 

The second derivative of the momentum dependent susceptibility at T = has a form 

flt^o, = 0) = 

x J d£p p {p: + £A) [(-3£ 4 + 5£ 2 - 8/3) exp (-3£ 2 /2) - Ei (-3£ 2 /2)] , (49) 

where = ±4t+8t' , Ei(a;) = dt e~*/i, p p {e) is the average value of square of electronic 
velocity in the isoenergetic surface + \i = e, 

X f dcj 

ft( £ ) = j\rE^ + ^" £k ^ = j 7^p Wk ( 50 ) 
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The calculation of the integral in fl50|) results in 



x 



it 2 (t 2 + at' - At' 2 ) 1 / 2 
-2\e\t + 4fe - e 2 )F(k 2 ) + 8 (V + et 7 - 4t' 2 ) E(/c 2 ) + 



x 



e 2 t 



where 



t + 2t'sign£: 
-At + e - 8t'){At + e - 8t') 



-n 1 



4t + 8t'signe ' 



16(t 2 + et'-4i') 2 

and F, E, II are the elliptic integrals first, second and third kind respectively. According to 
m> 9j#q=o(" = 0; A ^ 0) = -l/12 P ;'(/i). 

Plots of the functions of p p (e) and the second derivative of susceptibility with 
respect to momentum d 2 x -ff q= o (oj = 0; A — > 0) are presented in Fig. [131 The sign of 

(t/2^H q=0 (w=0;A^0) 
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Figure 13: The prot of : (a) p p (e);{b) d^Hq =0 {oj = 0; A -> 0) vs /x 

(9 2 ^iJ q=0 (ti; = 0; A — > 0) coincides with the sign of /x. Hence, the susceptibility of nonin- 
teracting electrons has its local maximum in the point q = for electronic concentration 
n < wvh and minimum for n > nyH- The second derivative of the susceptibility with respect 
to momentum tends to plus (minus) infinity as the electronic concentration approaches the 
Van Hove filling from the above (below). 
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